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1. Introduction 



The phenomenon of mixing in the r]-r]' system has attracted much attention since the 
advent of SUf{3) symmetry. Conventionally, one regards r/ and r/' as linear combinations 
of octet and singlet basis states, parametrized by a mixing angle Op. A determination of its 
value can be achieved either diagonalizing a mass matrix, as done in Chiral Perturbation 
Theory (xPT), or from phenomenology Several exhaustive analyses surveying many 
different processes have been performed along the years with values for the mixing angle 
ranging from —24° to —10° |2|-|10|. In all these previous analyses the dependence with 
energy of the r]-r]' mixing angle has been neglected. A recent analysis assuming such an 
energy dependence is also present in the literature |jll|. 

The r]-r]' mixing is strongly connected with the U{1) anomaly of QCD. Recently, it has 
pointed out that an extension of taking into account the effects of the chiral anomaly 
through a perturbative expansion leads to a description of the r]-rj' system in terms of two 



mixing angles [12|. In this framework, the four pseudoscalar decay constants associated 
to the matrix elements of axial-vector currents are written in terms of two basic decay 



constants /§, /o and two angles 9s, Oq |13|. Alternative approaches also including the chiral 
anomaly ab initio but using a more conventional one mixing angle scheme are discussed 
in Refs. P, 14, 15, 16, 17]. On the phenomenological side, the idea of using a two mixing 



angle scheme has been explored in Refs. [18, 19 1 



The aim of this work is to perform an updated phenomenological analysis of various 
decay processes using the two mixing angle description of the rj-rj' system. The analysis 
is carried out in the octet-singlet and quark-flavour basis and our theoretical predictions 
compared with the latest experimental data. The comparison will serve us to check the 
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validity of the two mixing angle scheme and its improvement over the standard one angle 
picture. The analysis also tests the sensitivity to the mixing angle schemes. 

In Section ^, we shortly introduce the notation used in the analysis. In Section |3|, we 
compute the radiative decays (??, V) ~^ 77 ^iid the ratio Rj/^ = ^{J/tp — > fl'l/vi) iii the 
two mixing angle scheme of both the octet-singlet and quark-flavour basis in order to obtain 
in each basis the preferred values for the mixing angles involved in the analysis. Section 
§ is devoted to the consequences of our approach in the context of the radiative decays 
of the lowest-lying vector and pseudoscalar mesons, V — > P7 and P V'j, respectively. 
In Section we compare our best results for the mixing parameters with the theoretical 
expectations of Large Nc xPT and previous phenomenological analyses. The implications 
of our analysis on these Large Nc predictions are also discussed. Finally, in Section 
|6|, we present our conclusions. 

2. Notation 

The decay constants of the rj-r]' system in the octet-singlet basis fp {a = 8,0; P = r],rj') 
are defined as^ 

{0\A';^\P{p)) = ^rpp,, (2.1) 



8 

where Af^ are the octet and singlet axial-vector currents whose divergences are 
df^A^^ = --^(m„m75U + mddi^^d - 2mssi'j5s) , 

where is the gluonic field-strength tensor and = igM'^"/^^'"^ its dual. The 



(2.2) 



divergence of the matrix elements (2.1) are then written as 

(OiaM^lP) = f},ml , (2.3) 

where mp is the mass of the pseudoscalar meson. 

Each of the two mesons P = r],r]' has both octet and singlet components, a = 8, 0. 
Consequently, Eq. ( |2.1] ) defines four independent decay constants. Following the convention 
of Refs. Ijl^, ^ the decay constants are parameterized in terms of two basic decay constants 
fs, fo and two angles 6*8, 6*0 



fv fv \ _ ( /s cos 6*8 -/osin(9o 
/°J ~ [fssmOs /o cos 00 , 



(2.4) 



This parametrization is the most general involving two independent axial-vector currents 
and two different physical states. 



^The axial- vector currents are defined as A"^ = 57^75^5 with the normalization convention fn 
V2F^ = 130.7 MeV. 
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(2.7) 



(2. 



Neglecting the contribution of the up and down quark masses, as in Ref. |14|, the 
matrix elements of the chiral anomaly between the vacuum and {r],r]') states are 

{0\^GG\v) = y|m2(/8cos08 - %/2/osin0o) , 

r- (2-5) 

{0\^GGW) = ^fm2,(/8sin08 + ^/2/ocoseo) • 

Analogously, in the quark-flavour basis the decay constants are parameterized in terms 
of fg,fs and (l)g,(j)s: 

( fv fv \ _ ( /acosc/)^ -/^sin(/>A ^^^^ 
V / V /g sin 0q fs cos 0^ / 

and the non-strange and strange axial-vector currents are defined as 

= = ^(^° - ^/24) . 

The divergences of these currents are 

d^'Al = ^/2(m„m75n + mdil^d) + ^^Gf.^G'"''^'^ , 

and therefore the matrix elements of the chiral anomaly in this basis are 

(OI^GGIt?) = A^2^^^os0, , (OI^GGIV) = -|m^,/,sin</-, . (2.9) 

3. Experimental values for the ^§,0 ^ind 0^ ^ mixing angles 

In order to reach some predictions from our two mixing angle analysis we must first know 
the values of 9% and preferred by the experimental data. We will use as constraints^ the 
experimental decay widths of {r],r]') ^ 77 Q 

r(f? 77) = (0.510 ± 0.026) keV , 

(3.1) 

r(??' ^ 77) = (4.29 ± 0.15) keV . 

Analogously to the vr'' — > 77 case, one assumes that the interpolating fields r/ and r]' can 
be related with the axial- vector currents (see e.g. Refs. |T^) in the following way: 

, , fO,d''Al{x)-f»,dt'AO(x) 

'/ l-^J ~ "Z;^ fO fS _ f8 fO • 



^We choose such constrains because those decays are well understood in terms of the electromagnetic 
anomaly (see e.g. Ref. pO|). 
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This leads to^ 

Y(r] 77) 
r(ry' ^ 77) 



967r^ 




1 a^mj^ 
1 ~ 967r^ 


V 






967r^ 


I fO f 8 fS fO 
y-lv-lr]' •'V-'r]'^ 


j ~ 967r^ 



(3.3) 



Because of the four unknown parameters {Og, Oq, fs and /o) that appear in Eq. ( |3.3| ), in order 
to get their allowed values we need two additional constraints (apart from the experimental 
constraints in Eq. ( p. ID ). On the one hand, we will use the well established prediction of 
xPT: /g = 1.28/7r (/tt = 130.7 MeV) as a theoretical constrain (see variants on this later). 
On the other hand, we will use the experimental value of the ratio 



R 



r(j/v; ^ V7) 



5.0 ±0.6 



(3.4) 



r( J/^ ^ r?7) 

According to Ref. [21|, the radiative J/tp P7 decays are dominated by non-perturbative 
gluonic matrix elements (see Ref. Q for further comments on the accuracy of this state- 
ment): 

{0\GG\r]') ^ 



R 



J/V 



(0|GG|r?) 



Pv 



(3.5) 



where pp = Mj/^{1 — mp/Mj^^)/2 is the three- momentum of the P- meson in the rest 
frame of the decaying J/il; (with mass Mj/^). Using Eq. (|2.5D one gets 



R 



m2,(/B,+v/2/0,) 



m2(/8+V2/0) 



(EnL 

V Pv 



ifs sin 6*8+^2/0 cos 60) 
m,2 (/8 cos 6*8 - V2fo sin 6*0 ) 



(3.6) 



Pv 



Comparing the experimental values of r(r/, rj' — > 77) and Rj/^p with the theoretical predic- 
tions shown in Eqs. (|3.3| ) and ( |3.6| ), one obtains 

^8 = (-22.2 ±1.8)° , ^o = (-8.7±2.1)° , /o = (1.18 ± 0.04)/, . (3.7) 



Using instead of fs = 1-28 fj^ the prediction of Large Nc xPT [g2|, g3|]: /s = 1.34/,^ ]13|], we 
would have found 



08 = (-22.9 ± 1. 



-6.9 ±2.0)°, /o = (1.20 ± 0.04)/, 



(3.8) 



The previous results are in agreement with the expectations of Large Nc xPT |jT2|, |T3[ 
and with phenomenological analyses [18, pi, |23] (see also Ref. pH] for a comparison of 



various analyses). Notice that the numerical value of /o cannot in principle be determined 
unambiguously, because it depends on the renormalization scale |2^. However, the 
scale dependence of the anomalous dimension of the singlet axial-vector current is very 



^Note that if one assumes a one mixing angle scheme {Os = 60 = 
e.g. Ref. @). 



the standard result is obtained (see 
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weak [12|. Therefore, in our numerical analysis we neglect that dependence and consider 
/o as a constant. The results of Eqs. ( p. 7] ) and ( p.^ ) constitute the first result of the present 
analysis. It is worth noting that the 6s and rnixing angle values are different at the 3a 
level, while the values of the pseudoscalar decay constant /o remain compatible with other 
results. 

In order to check the consistency of these results one can perform the same analysis 
but fixing 9% = 6q = 9, i.e. a one mixing angle scheme. Keeping /s = 1.28/7r and using 
^, V 11 in order to fit /o and 9, one gets: 9 = (-22.2 ± 1.7)° and /o = (1.07 ± 0.04)/^ 
—or 9 = (-22.9 ± 1.7)° and /o = (1.06 ± 0.04)/^ if /g = 1.34/^. In these cases, however, 
the predicted value for Rjj^ — 2.1 (or ~ 1.7) is incompatible with the experimental result 
in Eq. (0). Alternatively, leaving /§ as a free parameter one gets: 9 = (— 17.7±1.2)°, /§ = 
(1.02 ± 0.06)/^ and /o = (1.09 ± 0.04)/^. In this other case, Rj/^ ~ 3.5 and the value of 
/s is incompatible at the 3(T level with the predictions of xPT and Large Nc xPT. 

Therefore, a two mixing angle scheme leads to significantly different mixing angles 
in the octet-singlet basis to describe the radiative rj^Vj' — > 77 decays and the ratio Rjj^ 
simultaneously. 

In the quark-flavour basis (see the Appendix for the corresponding theoretical expres- 
sions), a simultaneous fit of r/, r/' — > 77 and Rj/^ with fq = f^^ suggests that, at the current 
experimental precision, (pq ~ 39.8° and (pg — 38.6°. Therefore, in this particular basis only 
one mixing angle seems needed to describe the data at the current experimental accuracy. 
If, for the exercise, we force this equality (which is not based in theory), (pg = cpg = (p, the 
result of the fit is 

/g = (1.07 ± 0.03)/^ , /, = (1.37 ± 0.27)/^ , (39.0 ± 1.7)° , (3.9) 



which is in agreement with earlier results [18|. The value of the mixing angle (p is fixed 



from the ratio Rj/^. Keeping however fq = /tt, as predicted by the OZI-rule, leads to a 
poor fit with f, = (1.39 ± 0.26)/^ and p = (39.8 ± 1.5)°. 

Let us mention at this stage that the fact of describing in terms of one or two angles is 
in itself not a relevant theoretical consideration, and that the resulting physical description 
should of course not depend upon the choice of the initial basis. In a general formulation, 
there are indeed two angles. Furthermore, the choice of the starting point (octet-singlet or 
magic mixing) is in no way obvious. The relevant symmetry, at least for the rj-r]' part of the 
present paper, is U{3), and subject to two types of breakings, one via the strong anomaly 
(inducing the S'C/(3)-singlet term, which can also be seen as an OZI-rule violation), the 
other via the strange quark current mass. The effect of both breakings is similar in size 
in the mass matrix, leading to a mixing structure somewhere between the two basis. The 
J/ip decays are of course also clear OZI violations in this framework. 



4. V-P electromagnetic form factors in the two mixing angle scheme 

In this Section, we want to extend our analysis to the V-P electromagnetic form factors. 
In particular, we are interested in the couplings of the radiative decays of lowest-lying 
vector mesons, V — > {rijri')^, and of the radiative decays rj' Vj, with V = p,uj,4>. In 
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order to predict such couplings we follow closely the method presented in Ref. Q where 
the description of the light vector meson decays is based on their relation with the AVV 
triangle anomaly, A and V being an axial-vector and a vector current respectively. The 
approach both includes SUf{3) breaking effects and fixes the vertex couplings gvP'y as 
explained below. 

In that framework, one starts considering the correlation function 

i I d'xe^'^^^Piqi + q2)\TJ^^ix)J^{0)\0) = e^,^pqU^^FvP,{qlql) , (4.1) 



where the currents are defined as 



J/?''^ = -^{uj^u =F dj^d) and = -57^5 



(4.2) 



The form factors values -Fyp^(0,0) are fixed by the AVV triangle anomaly (one V being 
an electromagnetic current), and are written in terms of the pseudoscalar decay constants 
and the (/)-uj mixing angle 9v as^ 

TP /n ^^ _ V3 ^ ■'v' 



V3 f^-V^ff, 
^ {cev-sev/V2)fO,-sdvff^, 



Fpri'y{0,0) — 4^ fO fH _ f8 fU ) 



7? (n - 1 {cev-sev/V2)f°-sevf^ ^ ' > 



p /-n m - 1 (^ev+cey/v/2)/0+eev/« 

Using their analytic properties, we can express these form factors by dispersion relations 
in the momentum of the vector current, which are then saturated with the lowest-lying 
resonances: 

Fyp^(0,0) = — WP7 + --- > (4.4) 

my 

where the dots stand for higher resonances and multiparticle contributions to the corre- 
lation function. In the following we assume vector meson dominance (VMD) and thus 
neglect these contributions (see Ref. 0] for further details). 



The (f)-Lu mixing angle in the octet-singlet basis is defined as 

(j> — cos Ovi^s — sin ^vt^o , 
uj = sin 6vi^a + cos Bvi^o , 

with Ljg = -^{uu + dd — 2ss) and uq = -^{uu + dd+ ss). Experimentally, the value of the vector mixing 

J — 11 

angle is measured to be above its ideal value and is 6v = (38.7 ± 0.2)° l29|. The value obtained from a 
quadratic Gell-Mann-Okubo mass formula for vectors mesons is 9v = (38.6 ± 0.4)° M. 
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The fv are the leptonic decay constants of the vector mesons, and defined by 

{0\J^\V{p,X))=mvfvei''\p) , (4.5) 

where my and A are the mass and the hehcity state of the vector meson. The fv can be 
determined from the experimental decay rates |jl] via 

T{V^e^e-) = ^a'^4, (4.6) 

6 my 

with cy = (-^, ^^) for y = yO, tj, </>. The experimental values are 

fpo = (221 ± 2) MeV , 

fuj = (180 ± 3) MeV , (4.7) 
= (239 ± 4) MeV . 

Finally, we introduce the vertex couplings gvP'y, which are just the on-shell V-P elec- 
tromagnetic form factors: 

{P{pp)\J^^\Vipv, A))|(p^_p^)2=o = W/3PpPy4(^) ■ (4-8) 

The decay widths of P — > V'j and V Pj are 

r(p ^ y^) = 1^2^^, 



7 \ mp 



^ Pj) = f^gl 



2 2 



mr, — m 



7 V "^l' 



3 



(4.9) 



V P gvPi (th.) gvp^ (cxp.) 



(1.59 ±0.11) GeV-i 
(1.35 ±0.06) GeV-i 
(0.46 ±0.02) GeV-i 
(0.46 ±0.03) GeV-i 
(-0.690 ±0.008) GeV-i 
(0.71 ± 0.04) GeV-i 

Table 1: Theoretical and experimental values of the on-shcU V-{ri,rj') electromagnetic vertex 
couplings in the octet-singlet 77-77' mixing angle scheme. For gvP'y (th.) we give the experimental 
errors coming from the decay constants /p,y and the mixing angle values 9s and ^o- We use 
9v — (38.7 ± 0.2)° for the cjj-u; mixing angle. Experimental values are taken from 0|. 

Eq. ([4.4D allows us to identify the gy couplings defined in (|4.8D with the form factors 
Fyp-y (0,0) listed in (|4.3|) . The couplings are expressed in terms of the octet and singlet 
mixing angles Os and the pseudoscalar decay constants /§ and /o, the cp-u mixing angle 



47r^/p cfloceg+aegsSo 

47r^/p cSoces+segseo 

m^^ (cev-s0v/s/2)c0o/fs-sevses/fo 

2V2TT^f^ cBocea+segsBa 

/^)39o//8+s9t.efl8//o 
2^7r2/^^ ceocSg+segsSo 

(39v+eev/y2)eeo//8+cev3e8//o 
2V2ir2/^ cSoceg+sSgseo 



= (1.48 ± 0.08) GeV-1 
= (1.23 ±0.08) GeV-1 
= (0.57 ±0.04) GeV-i 
= (0.56 ± 0.04) GeV-i 
= (-0.76 ±0.04) GeV~i 
= (0.86 ± 0.05) GeV-i 
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and the corresponding vector decay constants fy- These theoretical expressions are 
shown in Table [l|. We also include a numerical prediction for each coupling that should be 
compared with the experimental values extracted from (4.9) and Ref. [|[|. In the numerical 
analysis we have taken into account a value for the vector mixing angle of Oy = (38.7ib0.2)° 
p^ ]. Our predictions are obtained from the results in Eq. ( p.8[ ). The error quoted in Table 
|l] does not reflect the full theoretical uncertainty, but namely propagates the errors from 
( |3.8| ) and ( |4.7| ). The agreement between our theoretical predictions and the experimental 
values is quite remarkable with exceptions in the ujrj'j and 001]'^ cases. However, these 
two couplings merit some explanation. On one case, the experimental value for g^^^^^ has 
changed from (0.53 ± 0.05) GeV-^ of the PDG'02 [||] to the current (0.46 ± 0.02) GeV~^ 
due to the exclusion of the measurement based on e^e~ 777 by Dolinsky et. al. [pg|]. On 
the other case, the Qurj'-y coupling is rather sensitive to the (p-u: mixing angle; for instance 
setting 9y to the ideal mixing value of 35.3° reduces the coupling by a 10%. As seen from 
Table |l], the predictions for the decays (j) — > {r],r]')^, which are the best measured and are 
largely independent of By, are in good agreement with data. It is worth noting that the 
comparison shown in Table |l| is performed with /§ = 1.34/7^. Using instead /§ = 1.28/7r 
one obtains similar numbers with the exception of (7(/)»77 ~ 

(-0.80 ± 0.04) GeV-i and 

g(f>ri''Y = (0.91 lb 0.06) GeV~^ which are worse fitted. In this sense, the experimental data 
seem to prefer a value for /§ higher than the one predicted by standard xPT. Once again, 
the 2-angle calculation is the relevant one. For completeness however, we have included in 
the table below same comparison but fixing 9^ = 60 = 9- In this case, the results are very 
similar for /g = 1.28/,^ and /s = 1-34/7^ but in both cases the /d.ol. is increased by a 
factor of 3 as compared to the value in the two mixing angle scheme. In particular, g^ri'-y 
are fitted to (1.21 it 0.07) GeV~^ and (1.20 it 0.06) GeV~^ respectively, which are in clear 
contradiction with data. 

Table | constitutes one of the main results of our work. Our analysis shows that the 
assumption of saturating the form factors Fy by lowest-lying resonances is satisfactory 
(a conclusion already reached in Ref. Q), and that the rj-rj' system described in the two 
mixing angle scheme (octet-singlet basis) fits the data much better than the one mixing 
angle scheme does. To quantify this improvement, we have performed various fits to the 
full set of experimental data assuming, or not, the two mixing angle scheme of the r]-r]' 
system. The results are presented in Table |2|. To check the consistency of our approach we 
have extended the fit in Eq. ( p.7| ) or Eq. (|3.8D to include all experimental data that account 
not only for the decay widths {ri,r]') — > 77 and the ratio Rjj^ but also the radiative decay 
widths oi V ^ P7 and P Vj. The theoretical constraint /§ = 1.28/,^ or 1.34/^ is 
relaxed in order to test the dependence of the result on the value of this parameter. The 
experimental constrain 9y = (38.7 it 0.2)° is also relaxed to test the stability of the fit. 

As seen from Table |^, a significant improvement in the x^/d-o.f. is achieved when the 
constrain 9^ = 9q = 9 is relaxed (in the most favorable case the x^/d.o.f. is reduced by 
more than a factor of 3), allowing us to show explicitly the improvement of our analysis 
using the two mixing angle scheme with respect to the one using the one mixing angle 
scheme. In general, the fits with 1.34/,^ are slightly better than those with /g = 1.28/,^. 
To leave /g as a free parameter does not make any substantial difference in the one angle 
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Assumptions 




R6SllltS 


X /a.o.t. 


Assumptions 




Results 


X /a.o.t. 


and 60 free 


eg = 


(-22.5 ± 1.3)° 


42.3/6 


Bs = eo = B 


B = 


(-16.9 ± 1.2)° 


81.4/7 


fs = 1.28U 


Bo = 


- (-8.0 ± 1.4)° 




fs = 1.28/^ 


h = 


(1.15 ±0.03)/, 




{"^^ 7-1-0 9^° 


JO — 






(IR 7 -1- n 9^1° 

— yoo. 1 ^ yJ.^j 








and 60 free 


Ss = 


(-22.9 ± 1.3)° 


31.2/6 


6*8 = 6*0 = e 


B = 


(-16.7 ± 1.1)° 


77.8/7 


/8 = 1.34/^ 


Bo = 


= (-6.6 ± 1.4)° 




/8 = 1.34/, 


h = 


(1.16 ±0.03)/, 




By = (38.7 ±0.2)° 


fo = 


(1.23 ±0.03)/^ 




By = (38.7 ± 0.2)° 








and 60 free 


es = 


(-23.8 ± 1.4)° 


18.8/5 


Bs = eo = B 


B = 


(-16.5 ± 1.2)° 


77.4/6 


fs free 


Bo = 


= (-2.4 ± 1.9)° 




fs free 


h = 


(1.37 ±0.05)/, 




By = (38.7 ± 0.2)° 


h = 


(1.51 ± 0.05)/^ 




By = (38.7 ± 0.2)° 


fo = 


(1.17 ±0.04)/, 






fo = 


(1.29 ± 0.04)/^ 












Bs and 60 free 


Os = 


(-24.0 ± 1.6)° 


18.6/4 


Bs = eo = B 


B = 


(-15.7 ± 1.4)° 


76.1/5 


fs free 


Bo = 


= (-2.5 ± 1.9)° 




fs free 


fs = 


(1.37 ±0.05)/, 




By free 


fs = 


(1.51 ± 0.05)^ 




By free 


fo = 


(1.17 ±0.03)/, 






fo = 


(1.29 ± 0.04)/^ 






ey 


= (36.5 ± 1.8)° 






8v 


= (39.4 ± 2.2)° 













Table 2: Results for the 77-77' mixing angles and decay constants in the octet-singlet basis of the two 
mixing angle scheme (left) and in the one mixing angle scheme (right). For every fit, the theoretical 
assumptions taken, the set of numerical results, and the value of the x^/d.o.f. are shown in the first, 
second and third column respectively. The fitted experimental data includes the decay widths of 
(77, 77') — > 77, V P7, P V-f, and the ratio Rj/^. 



case while in the case of two mixing angles the value increases up to /§ = (1.51 it 0.05)/7r, 
in disagreement with the expectations of and Large Nc xPT. However, the analysis 
of the parameter correlation coefficients reveals a strong positive correlation between /§ 
and ^0 (+0.674) — the latter quantity being taken as a negative number in the present 
convention; in other terms, the correlation to the absolute value of is negative. The 
correlations between /§ and or and ^re negative and much weaker (—0.126 and 
—0.124, respectively). It is also observed that when /§ is fixed, the remaining correlation 
coefficients are smaller. This increase is translated into a considerably better fit. Finally, 
it is seen in Table ^ that relaxing the experimental constrain Oy = (38.7 it 0.2)° does not 
produce any effect on the fits. 

Up to now, we have shown in the octet-singlet basis the need for a two mixing angle 
scheme in order to describe experimental data in a better way. In the following, we proceed 
to perform the same kind of analysis but in the quark-flavour basis. We will see that at the 
current experimental accuracy, the two angles have compatible values in this case. There 
is however no strong reason to impose this equality as a constraint. See also above for a 
discussion of the respective effects of the strange quark mass and anomalies in the mixing. 

In Table |3| we present the results of various fits taking into account all experimental 
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Assuiiiptioiis 




RCSliltrS 


X /a. O.I. 


Assuiiiptions 


Results 


X /d.o.i. 


4)q and 0s free 


<Pq 


= (40.4 ± 1.2)° 


34.6/6 


<pq = (pa = <p 


<p= (40.8 ±0.9)° 


34.9/7 


fq ~ fn 


<Ps 


= (41.3 ± 1.3)° 




fq — /tt 


fs = (1.66 ±0.06)/^ 




(pv = (3.4 ±0.2)° 


fs = 


(1.66 ±0.06)/^ 




<pv = (3.4 ± 0.2)° 






<pq and (ps free 


<Pq 


= (39.9 ± 1.3)° 


18.8/5 


<pq = (ps = <p 


<p= (40.6 ±0.9)° 


19.4/6 


fq free 


4>s 


= (41.4 ± 1.4)° 




fq free 


fq = (1.10 ±0.03)/^ 




<t>v = (3.4 ±0.2)° 


fq = 


(1.09 ± 0.03)/^ 




4>v = (3.4 ± 0.2)° 


fs = (1.66 ±0.06)/^ 






fs = 


(1.66 ±0.06)/^ 










(pq and (ps free 


<Pq 


= (39.8 ± 1.3)° 


18.6/4 


<pq = (ps = <p 


<p= (40.4 ± 1.0)° 


19.3/5 


fq free 


<Ps 


= (41.2 ± 1.5)° 




fq free 


fq = (1.10 ±0.03)/^ 




<pY free 


fq = 


(1.09±0.03)/,r 




<pv free 


fs = (1.66 ±0.07)/^ 






fs = 


(1.67 ±0.07)/^ 






(pv = (4.1 ± 2.2)° 






<Pv 


= (4.2 ± 2.1)° 











Table 3: Results for the 77-77' mixing angles and decay constants in the quark-flavour basis of the 
two mixing angle scheme (left) and in the one mixing angle scheme (right). The conventions are 
the same as in Table 0. 



data available — namely {r],r]') 77, V P7, P — > Vj, and the ratio Rj/^ — and the 



corresponding theoretical expressions found in Eqs. ( A.2 )-( |Ay5| ) of the Appendix. The 
theoretical constrain fq = f^^ is relaxed to test the dependence of the result on the value of 
this parameter and the same happens with the experimental constrain (py = (3.4 ± 0.2)°. 
As seen from Table ^, there is no significant difference at the present experimental accuracy 
between the x^/d-o.f. of the fits when data are described in terms of two mixing angles (in 
the quark-flavour basis) or (j)q = (j)s = (p is imposed. It is however important to notice 
that the fit considerably improves when the parameter fq is left free. In this case the value 
obtained fq = (1.10 it 0.03) fj^ is incompatible with the value in the large Nc limit fq = f-,^. 



5. Discussion about the mixing parameters 

In this Section, we compare our best results for the pseudoscalar decay constants and mixing 
angles in the octet-singlet and quark-flavour basis with the theoretical expectations of Large 
Nc xPT and previous phenomenological analyses. The value of the mixing parameters 
extracted from our best fit in the two mixing angle scheme are 

/8 = (1.51 ± 0.05)/^ , ^8 = (-23.8 ± 1.4)° , 
/o = (1.29 ± 0.04)A , eo = (-2.4±1.9)° , 

in the octet-singlet basis, and 

fq = (1.09 ± 0.03)/^ , (t>q = (39.9 ± 1.3)° , 
fs = (1.66 ± 0.06)A , (Ps = (41.4 ± 1.4)° , 
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in the quark-flavour basis. These values are extracted from a comparison with experimental 
data only assuming that the pseudoscalar decay constants involved in the corresponding 
processes follow the parametrization given in Eqs. ( |2.4[ ) and ( |2.6| ), respectively. At the 
present accuracy, our results satisfy the approximate relations existing between the two 
different sets of mixing parameters^ |15]: 



/8 = ^1/3/2 + 2/3/, 
/o = J2/3/2 + 1/3/1 



2 

s ; 



arctan(\/2/s//g) , 
ar ctan{^/2 f g/ fs) . 



(5.3) 



Large Nc xPT predicts [12, |T 
J8 



f2 
JO 



/s/o sin (08 
in the octet-singlet basis and 



^+/^A, 



fn) 



(5.4) 



fqfs sin( 



s) 



3 •'T^'-i 



(5.5) 



in the quark- flavour basis |26|. These expressions are valid at next-to- leading order in 
the Large Nc expansion where the octet-singlet (and quark-flavour) pseudoscalar 

decay constants can be written in terms of the known /^r and fx decay constants and the 
unknown OZLrule violating parameter Ai. Using the experimental constrain fx = 1.22/7^, 
one obtains*^ fs = 1.28/^, = -20.5°, /o ~ L25, and 6^0 ~ -4°. Our results in 

Eq. ( |5.1| ) are quite in agreement with the former values except for the case of fg. Note 
however that in Ref. [12| the value of /§ is fixed from theory while in our analysis it is fitted 
from a direct comparison with experimental data where a positive correlation between /g 
and 9q appears. For /§ = 1.34/^ the results of the fit are in perfect agreement with the 
predictions from Large Nc even though the quality of the fit is slightly reduced. 

As seen from Table ^, if the constrain fs = 1.28/7r is imposed one gets a worse fit. 
The same kind of comparison can be performed in the quark-flavour basis. From the 

(L07 ± 0.02)/, 



phenomenological analysis |18] the values /, 



/, = (L34 ± 0.06)/, and 



(f) = (39.3 lb 1.0)° are obtained. In this case, our results in Eq. (5^) fairly agree with 
the exception of fs which clearly disagrees. This difference may be due to the fact that 
the analysis in Ref. ||l^ is based on a different set of experimental data not including for 
instance the very precise and recent (j) {r], rj')j decays which are very dependent on fs 
(see Eq. ([AlsD). 

Our best results for the mixing parameters in Eqs. ( |5.1| ) and ( |5.2| ) can be used to 
check the consistency of Eqs. ( |5.4D and (^.5|), and therefore to test the reliability on the 



The relations (5.3) are obtained once tlie non-strange and strange axial- vector currents are written in 



terms of the octet and singlet ones (see Eq. (2.7)) and are valid for 



A value of fs — 1.34/^ is obtained if chiral logs and higher order contributions are also taken into 
account lll3|. 
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Large Nc xPT framework. Accordingly, our fitted values for /g and /o together with the 
third equation in ( |5.4| ) can be used to get 9s — Oq = (—13.7 it 0.6)° as a prediction for 
the difference of the two mixing angles in the octet-singlet basis. If one compares this 
prediction with our result 6s — 0q = (—21.4 it 2.4)° a disagreement is again obtained. Using 
fs = 1.34/^ and our fitted value for /o one gets — Oq = (—16.2 it 0.4)° to compare with 
our result 9s — 9q = (— 16.3± 1.9)° (see Table ^. The second equation in (5^) may be used 
to get a prediction for the OZI-rule violating parameter Ai once /o is provided. Using our 
value for /o one gets Ai = 0.34 it 0.10. Ai can also be extracted from the set of equations 
dU). From the fitted values for fg^s and (pg^s the results Ai = 0.32 ± 0.10, 2.34 ± 0.60 and 
—0.10 lb 0.13 are obtained using the first, second and third equation in (|5.5D , respectively. 
Indeed, the same behaviour is also observed in the phenomenological analysis of Ref. |l8|| . 
The values obtained for the mixing parameters in this analysis lead to values for Ai that 
are incompatible with the constrain Ai = imposed in their fit. 

The numbers obtained in the former discussion hint at a disagreement with Large Nc 
xPT. Some care should be taken however to qualify this statement. On the one hand, 
the statistical significance is still limited, but will obviously improve as critical channels 
are measured more acurately. On the other hand, it is plainly clear that the critical 
information comes here from radiative decays of vector mesons, where we use vector meson 
dominance. Questions have been raised as to the consistency of vector resonances and 
VMD with Chiral Perturbation Theory and short distance QCD cf. e.g. |31]. This latter 
remark should however be moderated by the fact that our approach was carefully tested 
in the charged mesons sector (away from the eta mixing problems) [Q]. 



6. Conclusions 



In this work we have performed a phenomenological analysis on various decay processes 
using a two mixing angle scheme for the rj-r]' system. First we have used the radiative 
decays {r],r]') — > 77 together with the ratio Rj/^ to fit the values of the pseudoscalar 
decay constants and the mixing angles in both the octet-singlet and the quark-flavour 
basis. Using the description of vector meson decays in terms of their relation with the 
AVV triangle anomaly, a theoretical prediction for the gvP'y couplings have been derived. 
The agreement between our theoretical predictions and the experimental values is very 
good and can be considered as a consistency check of the whole approach. Second we have 
extend our analysis to include the V — > P7 and P — > decays in the fits showing that 
a two mixing angle description in the octet-singlet basis is definitely required in order to 
achieve good agreement with experimental data. On the contrary, in the quark-flavour 
basis and with the present experimental accuracy a one mixing angle description of the 
processes is still enough to reach agreement. Finally, we have compared our best fitted 
mixing parameters with the predictions from Large Nc xPT showing a possible discrepancy 
with this framework. Higher accuracy data and more refined theoretical analyses would 
contribute to clarify the preceding issue. 
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A. Theoretical expressions in the quark-flavour beisis 

In the quark-flavour basis, the interpolating fields 77 and rf are expressed as 



r]{x) 



fs f>l_f1 fs 



vi'{x) 



the {r],r]') — > 77 decay widths as 



1 ^d^AUx)-fpi^Al{x) 



(A.l) 



r(?? ^ 77) = ' — 



fs fQ_fQ fs 



Q^"^^ ( CqC4>s/fq-CsS4>Jfs'^ ^ 



rfa/ . - " c,r,-csf?, \ _^(c, 

with Cq = 5/3 and Cg = V^/S, and the ratio Rj/^ as 



32^ 
2^3 



CtpsCtp„ + S<PqS4 



qS4>s/ fq+CsC(i>q/ fs 
CCpsC<pq + S(pqS<ps 



(A.2) 



R 



(A.3) 



The form factors Fvp-y{0,0) are written in terms of the pseudoscalar decay constants 
(/?) fs) and the 4>-u mixing angle in the quark-flavour basis {(py = Oy — arctan -^): 

fs 



F<^^^(0,0) 



3 -^V 



1 c<t>vf^,-2s(f>vf!^, 



i^a,r,'7(0,0) = 

F<^^/^(0,0) = 



myfu;> 



Finally, the vertex couplings gy p-y in this basis are 



3m„ 



C4>s/fq 



3mp c(t> 



3mp S(j>s/ fq ^ 3mp s4> 

47r^/p c4>sC4>q+Scl>qS(t>s ~ iw^ fp fq 



C<t>VC<l>s/ fq-2s4>vs4>q/ fs 



iTp/uJ Ccj)sC4>q+s4>qS4>3 

- C(l>vS(j>s/fq + 2s(f>vc4)q/fs ^ 

in-'f^ C<j>sC<j>q + S<j>qS(f>s ~ 



— S<pvaj>s/fq+2c<pvS<j>q/fs 

y'I'Vl TPJ^ aPsC4>q+S<t>qS,Ps 

— S(j}vS<ps/ fq-2c<l>vC(f>q/ fs 

- TPJ^ C<j>sC<pq+S<PqS4>s 



^U<Py^ + 2s4>y'- 



S(t>Vj^ — 2c(l)y 



(A.4) 



(A.5) 
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where the approximations are valid for (j)q = (pg = <t>- 
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